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Abstract 

We propose and develop the complex scaled multiconfigurational spin-tensor electron propaga¬ 
tor (CMCSTEP) technique for theoretical determination of resonance parameters with electron- 
atom/molecule systems including open-shell and highly correlated atoms and molecules. The mul¬ 
ticonfigurational spin-tensor electron propagator method (MCSTEP) developed and implemented 
by Yeager his coworkers in real space gives very accurate and reliable ionization potentials and 
attachment energies. The CMCSTEP method uses a complex scaled multiconfigurational self- 
consistent field (CMCSCF) state as an initial state along with a dilated Hamiltonian where all of 
the electronic coordinates are scaled by a complex factor. CMCSCE was developed and applied 
successfully to resonance problems earlier. We apply the CMCSTEP method to get Be“ shape 
resonance parameters using lAsllpbd, 14sl4p2(i, and 14sl4p5(i basis sets with a 2s2p3(iCAS. The 
obtained value of the resonance parameters are compared to previous results. This is the first time 
CMCSTEP has been developed and used for a resonance problem. It will be among the most ac- 
cnrate and reliable techniques. Vertical ionization potentials and attachment energies in real space 
are typically within ±0.2 eV or better of excellent experiments and full configuration interaction 
calculations with a good basis set. We expect the same sort of agreement in complex space. 
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I. INTRODUCTION 


Resonances in electron-atom or -molecule scattering processes have attracted much at¬ 
tention. They play major roles in electron transport and energy exchange between electronic 
and nuclear motions, in vibrational excitation of molecules or molecular ions by electron im¬ 


pact, and dissociative attachmen 
damage by low-energy electrons 


:s and recombination [l, and as a mechanism for DNA 

,3. 


In order to avoid direct calculation of an outgoing wave in resonance problems, we use a 
complex coordinate scaling (CS) technique, which was proposed and developed by Aguilar, 
Balslev and Combes js, 6| and Simon j?] in the early 1970s. In this approach the electronic 
coordinates (r) of the Hamiltonian are scaled (or ’dilated’) by a complex parameter rj as 
r —)■ rjr, where r] = ae®® with a > 0 and 6 G (—tt, tt). Under this transformation, the bound 
states are real and are unchanged by complex scaling and the continua of the complex 
scaled Hamiltonian H is rotated by an angle 26 at each threshold such that the continuum 
states appear as complex eigenvalues of the complex scaled Hamiltonian H. The resonance 
parameters E = Ey. — hidden in the continua are exposed in complex space for some 
suitable rj, where Er and T^ are the resonance position and width of that resonance state, 
respectively. 

Other alternative methods have included the complex absorbing potential (CAP) BS 
instead of CS. CAP methods have not been shown conclusively to be superior to standard 
complex scaling. 

Previously, we developed the quadratically convergent complex scaled multiconhgura- 
tional self-consistent held [w, 11| (CMCSCF) method with step length control to obtain the 
resonance parameters. In real space, MCSCF with a small complete active space (CAS) has 
been proven to be a very effective method to describe nondynamical and some dynamical 
correlation correctly and is computationally cheaper than very large or full conhguration 
interaction (Cl) calculations 12| while still incorporating the fundamental physics of what 
is going on. Based on the CMCSCF initial state, we also developed a new method termed as 


the Ml method [ll|, ll3|, in which the comp 


of the M matrix dehned in MCSTEP 


14 


ex Ml matrix is constructed from the hrst block 


18| . This block allows for only simple electron 


removal and addition to orbitals with no more complicated processes allowed to mix in. 
MCSTEP, however, includes many additional operators which allow for more compli- 
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cated electron ionization and attachment processes to be included. MCSTEP is designed to 
calculate reliably the ionization potentials (IPs) and attachment energies (AEs) for atoms 
and molecules which cannot generally be handled accurately by perturbation methods. In 
addition to simple electron addition operators to all orbitals as in the Ml method, MC¬ 


STEP includes operators the allow for e 
to excited states within the CAS 


14 


ectron removal and electron addition to all orbitals 


18l |. In complex space, the Mi and CMCSTEP meth¬ 


ods use CMCSCF states as reference or initial state along with H. Both the CMCSCF 
and Ml methods have been previously efficiently used to study the ^PBe~ shape resonance 

3,HQ- 


Moreover, we have developed and implemented the complex scaled multiconhgurational 
time-dependent Hartree-Fock method (CMCTDHF) [also called the complex scaled mul¬ 
ticonhgurational linear response method (CMCLR)]. CMCTDHF uses CMCSCF state as 
the initial state. In real space multiconhgurational time-dependent Hartree-Fock (MCT- 


DHF) has 

properties lojj. CMCTDHF has previously been implemented and successully emp 
study Auger resonances for Li and Li-like cations [2^, and Be and Be-like cations 
Feshbach resonances for both Be'''(2p) 2^ and He(2s^) 23| systems, as well. 


jee n successfully used to study electronic excitation energies and linear response 

oyed to 
211 and 


In this work we implement the CMCSTEP method for the Be” shape resonance prob¬ 
lem using lAsllpbd, 14sl4p2d, and lAslApbd basis sets with a 2s2p3dCAS and compare 
our results with previous results. The reasons why we implement this method for resonance 
problem are (i) MCSTEP in real space works exceptionally well and gives the most accurate 
and reliable values of vertical IPs and AEs for general atomic and molecular systems, which 
are well consistent with experimental measurements 24j-l28|. so that we expect that CMC¬ 


STEP is able to give the most reliable values of resonance parameters; (ii) even though this 
approach has been implemented for the hrst time in complex space, this is a direct extension 
of the CMCSCF Q and Mi 0 methods which we previously developed and implemented, 
then we expect that the obtained results for this resonance problem will be different from 
those previously obtained and the most accurate HQ. 


The paper is organized as follows. In section H we discuss the theoretical part of CMC¬ 
STEP method. In section HI we present and discuss our results. Then conclusions follow. 
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II. THEORY 


The complex scaled electronic Hamiltonian is non-Hermitian. It is complex symmetric. 
This causes the wave function \'iprn) to be complex conjugate biorthogonal (CCBON) where 
= ^ij (* means complex conjugate) [29|. It is shown that creation operators are 



CMCSTEP IPs and AEs are obtained from the following the complex generalized eigen¬ 
value problem: 


MX/ = w/NX/, (1) 

where 

r 

x(2r+ l)‘''"(AfS„||{/!;(>),//,ftp( 7 p)}||AfS„), (2) 

and 

Nrp = 5^(-l)''“-^-''/-^’'lE(7.7p5o5o;P5/) 
r 

x(2r + l)'/2(Ars„||{/i;(7,),ftp(7p)}||AfS„), (3) 

ujf is and IP or AE from the At-electron initial tensor state |AtS'o) with spin 5*0 to At ± 1 
electron hnal ion tensor state | At± IS"/) which has spin Sf. IE is the usual Racah coefficient, 
hpi'jp) and h*p{^r) are tensor operator versions of members of the operator manifold with 
ranks 7 p and 7 ^, respectively, {,} is the anticommutator 

{A,B} = AB + BA, (4) 

and {,, } is the symmetric double anticommutator 

{A,B,C] = \{[A,\B,C] + {|A-S].C}}). (5) 
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CMCSTEP uses a CMCSCF initial state with a fairly small CAS and couples tensor ioniza¬ 
tion and attachment operators to a tensor initial state to a hnal state that has the correct 
spin and spatial symmetry even if the initial state is open shell and/or highly correlated. 


III. RESULTS AND DISCUSSION 


In this study, we investigate the low-lying ^PBe shape resonance problem using 


ACMCSCF, Ml and CMCS 


theoretically in the past 


33 


'FP methods. This resonance problem has been investigated 


43[|. Recently, we studied this resonance problem in terms of 
application of newly developed Mi method [ll]. The Be atom has a fairly large amount 


of nondynamical correlation because^ 
the principle ls^2s^ conhguration 


le conhguration has considerable mixing with 

151 1. so that both conhgurations need to be included 


nonperturbatively for accurate IP and AE calculations. 

It is a common practice to report the resonance energy relative to the total energy of the 
scattering target. In this work, in ACMCSCF calculations we report on the total energy of 
the continuum Be” species relative to that of Be atom as 


eAcMcscF(r?)=^f+'-<, (6) 

where, and Eq are total energies of the {N + 1) electron Be” resonance state under 

investigation and the N electron ground-state of the neutral Be atom, respectively, and 
subscript c and 0 refer to continuum and bound states, respectively. 

In ACMCSCF calculations we need to optimize each state separately, however, in Mi 
and CMCSTEP calculations we can obtain energies of all states simultaneously. In order to 
be consistent with the ACMCSCF calculation, we report on resonance parameter Ccmcstep 
obtained from CMCSTEP method: 


/ \ Q 

-CMCSTEP wJ ~ 


CMCSTEP 




N 

0 


( 7 ) 


where (^cmcstep ^ calculated from equation In the case of Mi calculations it is 

obtained from the Mi complex eigenvalue problem [ll| and we reported on results based on 
complex eigenvalues rather than i^cmcstep equation ([7]). 

For this resonance problem, Venkatnathan et al 33| found the 14sllp basis set to be the 
best one. We subsequently have shown that this 14sllp basis set is somewhat inadequate 
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for resonances and that at least 14p functions are much more reliable. Hence, we chose this 
basis set initially and added p and d functions to it using a geometric progression with a 
view to account for the diffuse nature of the resonances. Although for the very accurate IPs 
and excitation energies of Be, a larger 2s2p3s3p3d CAS which enables more correlation is 
necessary |^, we employ a 2s2p3dCAS with basis sets IdsllpSd, 14sl4p2(i and lAslAp^d in 
this calculation, since we have previously found that a larger CAS is unnecessary for accurate 


shape resonance calculations 


13l |. However, most of the IP basis sets are designed for IPs 


where tighter functions are necessary rather than for resonance calculations where what is 
needed are basis functions to describe the near continuum. First, we performed ACMCSCF 
calculations with all basis sets, and then followed it up with Mi and CMCSTEP methods. 
The hrst two methods have already been implemented for the resonance problem with other 
basis sets [w, 11|; however, here the CMCSTEP method is applied for the hrst time for this 
resonance problem. Of these, CMCSTEP should be the most accurate, efficient and reliable 
method. So far, there are no experimental results for resonance parameters of Be”. 

In Table I we present IP for the X ‘^S state of Be atom and AE for Be” resonance state 
obtained from ACMCSCF, Mi and MCSTEP (9 = Orad and a = 1) calculations, in which 
Im(E) = 0. A comparison of values of IPs and AEs to previously obtained theoretical values 
and experimental measurements presented in this table show that MCSTEP calculations 
are much better than MCSCF and Mi approaches, and will give more reliable values for 
resonance problems. 

In Table H we show a summary of the obtained values of Be” shape resonance for three 
different basis sets. In rows 2-4 of Table H we show results form ACMCSCF calculations. 
These give larger widths than Mi or CMCSTEP. Resonance parameters obtained from Mi 
and CMCSTEP methods shown in rows 5-10 are fairly consistent with each other, although 
CMCSTEP will be more accurate. The optimal values of a and 6 enables one to estimate 
the resonance parameters, and can be found by the system of equations below: 


dE _r]dE 
da a dp ’ 


( 8 ) 


dE dE 
-de = 

which form the trajectory method by determining E{aopt, 9opt) corresponding to the stability 
(loops, kinks, inflections, or any kind of ’’slow down”) in the plots of Im(E) as a function of 
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TABLE I. Ionization potential and attachment energies for X‘^S of Be atom and ^PBe ion, 
respectively. 


Method and basis set 

IP(eV) AE(eV) 

AMCSCF; 14sllp5d - 2s2p3d CAS 

8.478 

1.027 

AMCSCF: 14sl4p2d - 2s2p3ci CAS 

8.042 

1.124 

AMCSCF: 14sl4p5d - 2s2p3d CAS 

8.042 

1.121 

Mi: 14sllp5d — 2s2p3d CAS 

7.571 

0.836 

Mi: UsUp2d - 2s2p3d CAS 

7.576 

0.872 

Mi: 14sl4p5(i — 2s2p3d CAS 

7.576 

0.871 

MCSTEP: UsllpBd - 2s2p3dCAS 

9.508 

0.843 

MCSTEP: UsUp2d - 2s2p3dCAS 

9.506 

0.918 

MCSTEP: UsUpBd - 2s2p3dCAS 

9.506 

0.906 

Ref.[l^ -2s2pCAS 

9.50 


Ref. \15]—2s2v3s3v3d CAS 

9.31 


ExDt.f44l 

9.32 


1 a series of a {a trajectory) and a 

series 

of 9 values 1 


Vertical IPs and AEs in real space are typically within ±0.2 eV of excellent experiments 


and fnll confignration interaction calcnlations with a good basis set 


15j. We expect the 


same sort of agreement in complex space. Indeed this can be seen by looking at Table II 
where onr resnlts are in general agreement with considerably less accnrate methods. 

In Table III we have listed theoretical resnlts obtained by other workers. Onr cnrrent 
resnlts with onr best basis set (14sl4p5d) are quite far from complex Cl 39j and density 
functional theory (DFT) combined with a CAP ^ calculations. We note that the Cl calcu¬ 
lations did not include any effect of quadruple excitations. It is well known that these need 
to be included for accurate Cl energies and properties The complex DFT calculation 
contains parameters that are experimentally determined and also the basis set used is small 
(contracted Gaussian 5s4pld functions) and not adequate for resonance calculations. How¬ 


ever, the CMCSTEP obtained values in this work are 
by ASCF 38| and electron propagator methods 


airly comparable with those obtained 


33, 


m 


38, 


43l |. although these other 
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TABLE II. Summary of theoretical calculations for ^PBe shape resonance relative to 
ground state. 


Method and basis set 

(^opt 

Oopt(rad) 

Er (eV) 

F,. {eV) 

ACMCSCF: 14sllp5d-2s2p3(iCAS 

1 

0.49 

0.714 

1.541 

ACMCSCF: 14sl4p2d-2s2p3(iCAS 

1 

0.55 

0.816 

1.731 

ACMCSCF: 14sl4p5d-2s2p3(iCAS 

1 

0.55 

0.819 

1.736 

Mi: 14sllp5(i-2s2p3dCAS 

1 

0.36 

0.764 

0.796 

Mi: 14sl4p2(i-2s2p3dCAS 

1.03 

0.36 

0.790 

0.856 

Mi: 14sl4p5d - 2s2p3dCAS 

1.03 

0.37 

0.789 

0.874 

CMCSTEP: 14sllp5d-2s2p3dCAS 

1 

0.30 

0.768 

0.740 

CMCSTEP: 14sl4p2d-2s2p3dCAS 

1.03 

0.30 

0.795 

0.681 

CMCSTEP: 14sl4p5d-2s2p3dCAS 

1.03 

0.36 

0.756 

0.862 


methods will not be as accurate, since they are based on a single conhguration and Be atom 
is inherently multiconhgurational with the ls^2p^ conhguration mixing in strongly (i.e. 10%) 
with the ls^2s^ conhguration for the initial state. From our previous experience with calcu¬ 
lations for IPs and AEs for atomic systems 25|, we know that MCSTEP works very 

well, therefore we can say that value of resonance parameters obtained from CMCSTEP in 
this work are reliable. 


In Figure 1 we show the ^-trajectories for ^PBe~ shape resonance obtained from the 
CMCSTEP method. Curves shown in panels (a), (b) and (c) are corresponded to calculations 
with basis sets Idsllphd, 14sl4p2d and 14sl4p5(i with 2s2p3(iCAS, respectively. Crosses 
on each trajectory show a stabilized point. All trajectories show resonance points clearly 
along with an increased density of points. In all trajectories 9 starts a.t 6 = 0.01 rad at the 
top and increases to down with a step of 0.01 rad. 


Although we have here presented results for resonance parameter for an atomic system, 
^FBe~, the method can be implemented for investigating shape resonance parameters for 
molecular systems. We had shown for the ^If^N^ shape resonance using Mi, and it 
is quite consistent with previous literature results 45l-l48l| and experimental measurements 
49, ^]. In the molecular case 1^, the CS technique for the electron-nuclear Coulomb 
interaction potential — 2’/|r — R| has been implemented so that —{Zri~^)/\r — where 

Z is a nuclear charge, and r and R are the electronic and nuclear positions relative to an 


origin of a hxed molecular coordinate system 


341. We will report CMCSTEP calculations 
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TABLE III. Theoretical calculations for ^PBe shape resonance. 


Method 

Er (eV) 

F, (eV) 

Static exchange phase shift [3^ 

0.77 

1.61 

Static exchange phase shift plus polarizability phase shift [3^ 

0.20 

0.28 

Static exchange cross-section [3^ 

1.20 

2.60 

Static exchange plus polarizability cross-section [3^ 

0.16 

0.14 

ASCF with complex 14sl6p Gaussian basis set ]37} 

0.70 

0.51 

ASCF with complex 5sllp (Slater-type) basis set [38] 

0.76 

1.11 

Single, doubles, and triples complex Cl [39] 

0.32 

0.30 

S matrix pole {Xa ) [^, 

0.10 

0.15 

Complex density functional theory 

0.580 

0.223 

Second-order dilated electron propagator based on real SCF 

0.57 

0.99 

Bi-orthogonal dilated electron propagator (bases set 14slip ) [^: 



Zeroth order 

0.62 

1.00 

Quasiparticle second order 

0.61 

1.00 

Second order 

0.48 

0.82 

Quasiparticle third order 

0.54 

0.82 

OVGF third order 

0.54 

0.78 

Third order 

0.53 

0.85 

ACMCSCF (14sllp-2s2p3s3pCAS) [10] 

0.73 

1.58 

Ml (14sllp3(i-2s2p3s3p3dCAS) [Uj 

0.72 

1.12 

This work: 



ACMCSCF (14sl4p5d - 2s2p3dCAS) 

0.819 

1.736 

Ml (14sl4p5(i-2s2p3(iCAS) 

0.789 

0.874 

CMCSTEP (14sl4p5d-2s2p3(iCAS) 

0.756 

0.862 


for several molecules using this procedure in the near future. 
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(a) (b) (c) 

FIG. 1. The 0-trajectories for Be” shape resonance obtained from the CMCSTEP method. The 
curves shown in panels (a), (b) and (c) correspond to basis sets lAsllpbd, 14sl4p2(i and lAslApbd 
with 2s2p3(iCAS, respectively, and a cross shows a stabilized point. Computational parameters 
are a = 1 (a), 1.03 (b,c) and A0 = O.Olrad. 


IV. CONCLUSIONS 


In this work we have developed the CMCSTEP method and presented theoretical calcu¬ 
lations for Be” shape resonance using three different (ACMCSCF, Mi and CMCSTEP) 
methods. In our group we previously developed ACMCSCF and Mi methods, however, 
we here have developed and implemented for the hrst time the CMCSTEP method for res¬ 
onance problems, using three different bases sets IdsllpSd, 14sl4p2(i and 14sl4p5(i with 
2s2p3(iCAS. In CMCSTEP calculations we use the CMCSCF state as an initial state. The 
obtained values of ^PBe” shape resonance from CMCSTEP method are compared with 
previously obtained results in the literature. Based on our previous results from MCSTEP 
calculations for IPs and AEs for atomic and molecular systems, the results from CMCSTEP 
calculations are probably the most reliable and practical for resonance problems. 

In next step of our research work we intend to apply the CMCSTEP method to resonance 
problems for open-shell atomic and molecular systems, as well as studies for Feshbach and 
Auger resonances. 
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